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Abstract. The paper focuses on discrete-type approximations of solutions to 
non- homogeneous stochastic differential equations (SDEs) involving fractional 
Brownian motion (fBm). We prove that the rate of convergence for Euler 
approximations of solutions of pathwise SDEs driven by fBm with Hurst index 
H > 1/2 can be estimated by 0(S 2H ~ 1 ) (S is the diameter of partition). For 
discrete-time approximations of Skorohod-type quasilinear equation driven by 
fBm we prove that the rate of convergence is 0(S H ). 



Many equations which arise in modeling of processes in physics, chemistry, bi- 
ology finance, contain randomness. This randomness is not always well modeled 
by the classical Gaussian white noise (Brownian motion) because of long-range de- 
pendence, or long memory, of the processes under consideration. In this case the 
appropriate model for the randomness is fractional Brownian motion. Recall that 
B = (B t )t>o is called fractional Brownian motion on a complete probability space 
(f2, P) with Hurst parameter H £ (|, 1) if B is a centered Gaussian process with 
stationary increments and covariance Rh (t, s) — E(BtB s ) — \{t 2H +s 2H — \t— s\ 2H ). 

Numerical solution via time discretization of SDEs driven by Brownian motion 
has long history. We refer to the monograph |6J, which contains almost complete 
theory of numerical solution of such SDEs with regular coefficients. The paper 
[7j is devoted to Euler approximations for SDEs driven by semimartingales. Con- 
cerning numerical solution of SDEs driven by fBm, we mention first the paper [3], 
where equations with modified fBm that represents a special semimartingale are 
studied (recall that fBm itself is not a semimartingale) . Papers [TUl [H] study Euler 
approximations for homogeneous one-dimensional SDEs with bounded coefficients 
having bounded derivatives up to third order, driven by fBm, and prove that error 
of approximation is a.s. equivalent to <5 2ff-1 £t, and the process £t is given explicitly. 
These papers also discuss Crank-Nicholson and Milstein schemes for SDEs driven 
by fBm. The methods used by the authors cannot be applied to our case, because 
they require homogeneity and high regularity of the coefficients. 

We consider the stochastic differential equation on R d 



(1) X l t =X^ + J2 <r ij (s,X s )dBj+ Vis.Xjda, i = l,...,d, t € [0,T] 
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where the processes B l ,i — l,...,m are fractional Brownian motions with Hurst 
parameter H, Xq is a d-dimensional random variable, the coefficients a l \b l : £1 x 
[0, T] x R d — > R are measurable functions. 

The integral in the right-hand side of Q can be understood in the pathwise 
sense defined in [T2J, [11] or in Wick-Skorohod sense pQ . We treat the pathwise case 
first. We remind that the pathwise integral w.r.t. a one-dimensional fBm B can be 
defined as 

[ fdB = f {D a a+ f){s){Dlz a B b J){s)ds, 



where 



and 



(D« + f)(s) = 



1 
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are fractional derivatives of corresponding orders, 

B 6 -(«) = (B.-B 6 )l(.,6) 



(s). 



The integral exists for any a G (1—H, v) if, for example, / € C v (a, b) with i^+iJ > 1. 
Moreover, in this case pathwise integral admits an estimate 
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where Cq(u>) = C 


' SU Pa<s<fc 


\ D 



ds 



\m-f(u)\ 

(s - u) a+1 



du ds 



Denote a = (a t:l )dxm,b = (& l )dxi and for a matrix A — (a lJ )dxrn, and a vector 



y = (y*) dxl denote \A\ = £. j \a»\, \y\ = 

We suppose that the coefficients satisfy the following assumptions 

(A) a(t, x) is differentiable in x and there exist such M > 0, 1 — H < (3 < 



1, 



H 



!<«;<! and for any A > there exists such Mjy > that 



(B) 



1) \a{t,x)-a(t,y)\ < M\x - y\, x,y G R d ,t G [0,T]; 

2) |a xi o-(t,a;)-a !e4 <T(t,tf)|<M JV |a:-y|Ma;|,|i/| < A, ie [0,T]; 

3) |o-(t, x) - u(s,a;)| + \d Xi a(t,x) -d Xi a(s,x)\ < M\t-sf,xe R d ,t,s G 
[0,T]. 

1) for any TV > there exists Ln > such that 



|&(t,a;)-&(t,y)| <L N \x-y\, \x\, \y\ < N, t G [0, T]; 

2) |6(t,z)| <Z(l + |z|). 

As it was stated in [IT], under conditions (A)-(B) the equation ([T]) has the unique 
solution {X t , t G [0, T]}, and for a.a. u60 this solution belongs to C H_P [0, T] for 
any < p < H. Now, let t G [0,T],<5 = r„ = ^ = nS, n = 0, A. Consider 
discrete Euler approximations of solution of equation ([I]) , 



Z T n +l J T, 



+ 6 l (r„, 1^)5 + ^ a« (r„, ^JAB^ 
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and corresponding continuous interpolations 

m 

(3) Y t l > 5 = Y% + ft*( Tn , Y T s J(t - r n ) + ^{r n , Y t 5 J(B{ - <), t € [r„, r n+1 }. 

3 = 1 

Continuous interpolations satisfy the equation 

/•t m »t 

(4) y;- 5 = X* + b\t u , Yl)du + / * ij (tu, YfjdBl 

Jo - =1 Jo 

where <„ = t„ u , n n = max{?i : r„ < u}. 

For simplicity we denote the vector of solutions as X t = (^Q)i=i,...,d, vector 
of continuous approximations as Y t s — (y/'*)i=i,. ..,<*. Throughout the paper, C 
denotes a generic constant, whose value is not important and may change from line 
to line, and we write C(-), if the dependence on some parameters is crucial. 

The paper is organized as follows. Sections 2 and 3 are devoted to equations 
with pathwise integral. Section 2 describes the growth and Holder properties of 
approximations Y t s . We use here growth and Holder estimates of solution of corre- 
sponding pathwise equations from the paper [TT]. Section 3 contains estimates of 
rate of convergence for Euler approximations of the solutions of pathwise equations. 
It is well-known that in the case H = |, when we have SDE with Ito integral with 
respect to Wiener process, the rate of convergence is 0{5 1 / 2 ). It is natural to expect 
that in our case the rate of convergence might be 0(5 H ). Nevertheless, the rate 
of convergence is only of order 0(S 2H ^ 1 ) unless the diffusion coefficient is constant 
(where the rate is 0(S H ), as expected). Section 4 is devoted to discrete-time ap- 
proximations of solutions of SDEs with stochastic divergence integral with respect 
to the fBm (otherwise known as Skorohod integral, or fractional white noise inte- 
gral). It is shown that the rate of convergence is 0(S H ) in this case. The better rate 
of convergence is mainly because in the pathwise case there is no "Ito compensator" 
for the integral. 

2. Some properties of Euler approximations for solutions of 

pathwise equations 

In this section we consider growth and Holder properties of the approxima- 
tion process {Y t s ,t € [0,T]}. We need some additional notations. Denote <p u<v := 
\Y t S u -Y v s \ (« - v)-"- 1 for < v < t u < T, < a < 1, X* := sup < s < t '\X S \, 

Yf'* := sup 0<s<t |y/|. Further, for any < p < H there exists such C = C(ui,p) 
that for any < v < u 

(5) \B U -B V \ <C(^,p)( U - W ) H ^. 
We shall use the following statement [Til Lemma 7.6] 

Proposition 1. Let < a < 1, a, b > 0, x : R + — > R + be a continuous function 
such that for each t 

x t <a + bt a I (t- s)- a s- a x s ds. 
Jo 

Then x t < ac a e^{d a tb l l^~ a ^}, where c a = -ie 2 -^^-, d a = 2(T(l - a)) l/{1 ~ a) , 
r(-) is Euler 's Gamma function. 

We also establish technical lemma, which will be used later. 
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Lemma 1. There exists such C — C a > that for any s € [0, T], s ^ t s and 5 < 1, 
a E (0, 1) it holds 



J := ['(s-u)- 01 - 1 ( " {v-t v )- a dvdu<C5- a . 

JO Ju 

Proof. Evidently, 

J= (v-t v )- a (s-uy^dudv <a- x (v -t v )- a {s-v)- a dv 

JO JO JO 

Let i, = nS for some < n < N. Then 



rt 3 n-2 pT k + 1 A2n-l)8/2 rnS 

/ (v-t v )- a (s-v)-<*dv = Y, + + 

Jo k=0 Jrk J(n-1)5 J(2n-l)5/2 



We estimate the integrals individually: 



rl-2a 



f"k+l f T k + l 

/ <{s-T k+1 )- a (v-t^dv^il-ay^s-n+^S 1 - 

Jr k Jr k 
r(2n-l)8/2 /-(2n-l)5/2 

/ <(S/2)- a (v — t v )~ a dv < CS 1 

J(n-1)6 J{n-l)S 
rn8 rn8 

/ < {S/2)- a / (s-v)- a dv < CS 1 

J(2n-l)8/2 J(2n-\)8I2 



Therefore 



"-2 „ n 8 
j < CS l-2a + § -a yV _ Tk+i y*S < CS 1 ' 2 " + S~ a / (S - v)- a dv 
fc=n JO 



fc=0 

< C(5 1_2a + CV5~ a < CS- a . 



□ 
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Theorem 1. (i) Let the conditions (A)-(B) hold and 
(C)l) \*(t,x)\<C(l + \x\). 

Then for any e > and < p < H there exists 5q > and Q e ^s ,p C £1 such 
that P(Cl ei s 0tP ) > 1 — e and for any uj £ Q St s , p , 5 < 5 one has \Y t s \ < C{uj), 
\Yl - Y t 5 T \ < C(u)(t s - U) h -p, 0<r<s<T. 

(ii) If, instead of (A), 2) and (C) we assume that b and a are bounded functions, 
then |F/| < C(lu), \Y s s -Y r s \< C(lu){s - r) R -P, < r < s < T. 

In both cases C(uj) does not depend on 8. 

Proof. We can assume that S < 1. It follows immediately from (A), 1) and 3) and 

(0 that for any a G (1 - H, (3 A 1/2) 

(6) 









+/' 




Jo 




w 


Jo 



3=1 " U 

+ jdu + CoHV / \a l \t Ul Y t s )\u- a du 

m „t „ r 

J=1 Jo Jo 

|A*| + (C (u) + LT) + (C*oH + CT a ) \Y£ I u~ a du 
+ MC (u) J J* ({t r - t u f + \Y£ -Yf\ + |Y,f - Y t s u \){r-u) 



du dr. 



(We use here the equality t r = t u for t r < u < r.) Denote Ci(u>) := m(Co(u>) \_ a + 
LT) + \X \, C 2 (uj) := m(C (uj) + CT a ). Further, note that t r - t u < r -u + 6. 
Also, it follows from representations J3J that for any p G (0, H) 



(7) 



|y„ 5 - y£ | < l(i + \Y t 5 a \)(u-t u ) + c- c(w, P )(i + \Y t 5 u \)(u- t u ) 

<c 3 (Lu)(i+\Yi\y u -t u ) 



H-p 



\H-p 



where C 3 (w) = LT l ~ H ~ p + C ■ C{u, p). 
Moreover, for f3 > a 

P t := [ f "(U - t u f(r - uY^dudr < f I ' ' ((r - u)" + <5 /3 )(r - u)""" 1 ^* 

Jo Jo Jo Jo 

^{p-a)- 1 f ^-"dr + a^S 13 f (r-t r y a dr, 



and for any k > and any power n > — 1 



(r - t r ydr = (r - Tfefdr = CiS n+1 with Ci = (tt + 1)" 

whence 

(8) / (r-t r )" Q dr< / (r-t r y a dr = dNS 1 ' 01 = Ci<T a . 
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Therefore 

(9) P t < dT^ a+1 + ar l GiS f> ~ a < dT p - a+1 + oT x C x =: C 2 . 

Estimate now 



Qi 



\Y°-Y° (r -u)- a - L dudr, 



using and ©: 



(10) 



<C 3 (u)(l+Y t d '*)6 



S '*\*H-p-l 



(r - t r y a dr < Ci(uj){l + Y t d '*)S 



8,*\rH-a-p 



withC 4 (w) = C z {u))a- X -Ci. Note that Y/'* := sup < s < 4 \Y S S \ < oo for any* e [0,T] 
a.s. Substituting ([9]) and (fT0|) into ((6|), we obtain that 



(11) 



F/| <c 5 (w) + c 2 (w) / |y/J u - Q du + c 4 M(i + y/'*)(5 



+ C 6 (oj) / / tp r ^ u dudr 
Jo Ja 



r ^:*\ zH—a—L 



with C^{uj) = Cs(oj) + MCo(uj)C2, Cq(uj) = MCo(w). To simplify the notations, 
in what follows we remove subscripts from C(to), writing C(uj) for all constants 
depending on u>. 
So we can write 



(12) y/>* < C(w) (l + Y t s '*S B - a - p + J* \Yl 



u a du + 



ipr yU du dr 



o Jo 



In turn, we can estimate J Q S ip s u du. At first, similarly to the previous estimates, 
(13) 



K - Y '\< [ / " (l + \Yl\)dv+ / 3 (l+ \Yt \){v - u)- a 



*dv 



t 3 r t v 



\a(t v ,Y t s J - a(t z ,Y t s j\ (v - zY^dzdv 



< C(lo) \(t s - u) 1 -" + / |F/ I (v - u)- a dv + S (v- t v y a dv 



ip v ^ z dz dv + 



u J u 



\Y Z S -Y t s \(v- z)- a - x dzdv 



U J U 



multiplying by (s — u) a 1 and integrating over [0, t s ], we obtain that 



(14) 



pt 3 5 

/ f s , u du < c(uj) qi, 

J ° i=l 
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where 

(15) Ql 

(16) Q 2 S 



{t s - U y~ a (s~u)- a - i du< / {s-u)- 2a du < C; 



- «)~ 



\Y t s v | (v - u)- a dv 



(v - u)- a (s - uy^dudv < C / \Y t s I (s - v)- 2a dv 



where Co = f^°(l + y) a 1 y a dy; according to Lemma [T] 

Q 3 S :=S I '' '(s-u)- 01 - 1 I ' ' (v-t v )- a dvdu 
(17) Jo Ju 

< C5 f) 5- a < C. 

Further, using estimates ([7]), we can conclude that 
(18) 



< 



(s - u)~ 

t 3 pt v pzAv 



ip ViZ dz dv du 



o Jo Jo 



Vv,z{s — u) a 1 dudzdv<C I (s — v) a (p v ^ z dzdv. 



At last, using estimates (7) and Lemma [TJ we can conclude that. 
(19) 



(s-u)- 



U J U 



\Y* -Y t s Uv-z) 



-a-l 



dz dv du 



t 3 



<C(u)j (s-u)-*- 1 J j (v-z)- a - 1 dzdvdu-5 H - p (l+ Y t S J* } 
< C{u)(l + 



§H-p-ct 



Now, denote ip s :~ Y^'* + J Q S ip s ^ u du. Note that the integrals Q l s are finite for 
s = kS, i.e. for any s € [0, T], including s = t s . Then it follows from (fT2|) and 
(O-dlll) that 

^ t < C{uj) (l + Y t 5 '*8 H - a -P + J ((* - v)~ 2a + v- a )^ v dv) . 

Let e > be fixed. Note that all constants C(u>) are finite a.s. and independent of 
5. Thus, we can choose <5 > and ^l e ,s ,p sucn that C(ui)5 Q H ~ a ~ p < 1/2 on ^l e ,s ,p 
and P(fl et s , P ) > 1 — £■ Then for any w G £l e ,s a ,p 



whence 



i) t < C(u) + -fa + C(w) / {{t - v)- 2a + v- a )i> v dv, 



ipt < C(u) (l + t 2a j (t- v)- 2a v- 2a i) v dv^ , 



and it follows immediately from the last equation and Proposition!]] that ip t < 
whence, in particular, \Y t s \ < C(u), t £ [0, T], and J Q S ip u du < C(u). Moreover, 
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from (fl3|) with u = t r , r < s, taking into account that J t * s («— t v ) a dv < 8 a (t s —t r ), 
we obtain the estimate 



I - F/J < C(u;)((t s - t r y- a + 5P- a (t s - t r ) + (t s - t r ) 

+ 5 H -p f\v- t v y a dv} < c(«)(t, - t r y- a , 



and the statement (i) is proved, (ii) Let \b(t,x)\ < b, \a(t,x)\ < a. Then it is very 
easy to see that the estimate will take a form 



\Y t s \ < C(u)(\ + I j <Pr,ududS), 



13| will perform to 



;*|<c( w )((t s - ) i) 1 - a + (A ,3 + ^) f\v-t v y a dv 



(p ViZ dz dv 

and instead of (JT3I) — (JTHJ) we obtain 

tp s ,udu < C(lj)(i + J (s - v)~ a J tfv^dzdvj, 

whence the proof easily follows. □ 

3. The estimates of rate of convergence for Euler approximations 
of the solutions of pathwise equations 

Now we establish the estimates of the rate of convergence of our approximations 
(HJ) for the solution of the equation ([1]) with pathwise integral w.r.t. fBm. We 
establish even more, namely, the estimate of convergence rate for the norm of the 
difference X t — Y t s in some Besov space, similarly to the result of Theorem 1. Denote 

A UiS (X,Y 5 ) := \X S - Y s s -X u + Y*\ 

and assume for technical simplicity that = L, Mm = M in (A) and (B). 

Theorem 2. Let the conditions (A)-(C) hold and also 

(D) 1) Holder continuity of the coefficient b in time: \b(t,x) — b(s,x)\ < C\t — 
a\i, C > 0, 2H- 1< 7 < 1; 
2) the exponent [3 from (A) 3) satisfies (3 > H . 

Then: 

(i) for any e > and any p > sufficiently small there exists 8q > and Q, £t s . P 
such that P(Q e j ,p) > 1 — £ and for any u € Sl e ,& ,p> $ < $o 

U S := sup (\X S -Y S 5 \+ f' \A u . s (X,Y s )\(s-u)- a - 1 du) < C(lj) ■ S 211 ' 1 -", 

0<s<T V Jo 1 

where C{u>) does not depend on 8 and e (but depends on p); 

(ii) if in addition, the coefficients b and a are bounded, then for any p G (0, 2H — 
1) there exists C(u>) < oo a.s. such that Us < C{uj)8 2H ~ 1 ~ p , C{lo) does not depend 
on 8. 
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Proof, (i) Denote Zf := sup < s < f \X a -Y a s \. Then 
(20) 

Zf:= sup |X,-y/|< sup / \b(u,X u )-b(t u ,Y t s J\du 

0<s<t 0<s<tJo 

+ sup V I / (a^(u,X u )-a»(t v ,Y t 5 J)dBi\< / \b(u, X u ) - b(u, Y*)\du 

0<s<t ij=1 Jo Jo 

+ [ t \b(u,Y«)-b(t u ,Y l s l )\du+ [ \b(t u ,Y*) - b(t u ,Y t s J\du 
Jo Jo 

+ sup J2 I / (°- l3 (u,X u )-o-^(u,Yt))dBi\ 

0<s<t . Jo 
- - i,3 = l 

m „ s 

+ sup J2 I / (° l Hu,Y£)-v lJ (tu,Y*))dBl\ 

0<s<t- ■ , JO 
L >J — J- 

m „ s 6 

+ sup i / (^'(^^-^(^^Pii-E^ 

Now we estimate separately all these terms. Evidently, 

(21) h<L [ Z & u du. 

Jo 

Condition (D) 1) implies that for 5 < 1 

(22) 1 2 <C f \u- t u y du < C<S 7 < C5 2H -\ 

Jo 

As it follow from Theorem 2.2, for any e > and any p S (0, H) there exists So > 
and r^A.p C £1 such that P(ri e: 5 0-(9 ) > 1 — e and C(lu) independent of e and 5 such 
that for any ui 6 ^e,i5o,p ^ holds \Y t s —Y s s \ < C(uj) \t — s\ ~ p . In what follows we 
assume that 6 < Sq < 1. Therefore 

(23) h < L ■ C{uj)5 R - p ■ t < C(uj)6 H - p , uj G fie,* .P- 
Now we go on with I4. It follows from @ that for 1 — H < a < 1/2 



m r r 
h < C(u) Y / 



m r ft 

\<J ij (u, X u ) — a 1 3 (u,Y t 5 )\u- a du 



( 24 ) + / r \a^(r,X r )-a^(u,X u )~a^(r,Y r s ) + a^(u,Y^)\ 

Jo Jo 

x (r - uy a ^dudr =: I 7 + 1$. 

Evidently, 

(25) h < C(w) / ZfM- Q dw. 

Jo 
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According to (TTJ Lemma 7.1], under condition (A) 

|er(ii, xi) - a(t 2 ,x 2 ) - a(h,x 3 ) + er(f 2) Xa)\ < M \xi -x 2 - x 3 + x±\ 

(26) 



M \xi - x 3 \ ( \t 2 - tif + \xi - x 2 \ K + \x 3 - x 4 \ 



Therefore, Is < Ylk=s^t wnere 



t PT 

h = C{u) I / \X r -Y r s \ {r-uf- a - l dudr, 



IQ JO 
ft rr 



I m = C(uj) / \X r -Y r d \ \X r -X u \ K {r-uY^dudr, 
Jo Jo 



t rr 



Iu = C{u) / \X r -Y°\ \Y°-Y£\ (r-u)- a - L dudr, 
Jo Jo 

h 2 = C(w) [ f A^iX^Y^ir-uy^dudr. 
Jo Jo 

Taking into account that (3 > H > a, we obtain that 

(27) h < C(uj) [ Z & u du. 

Jo 

As it follows from 11, Theorem 2.1], under assumptions (A) and (B) for any < 
p < H there exists such constant C(u) that 

(28) sup \X t \<C{uj), sup \X t -X„\ <C{u)\t- s\ H ~ p . 

0<t<T 0<s<t<T 

Moreover, we can choose p > and a > 1 — H such that n{H — p) > a and 
H — p > 2H — 1, because nH > 1 — H. In this case 

ft rr pT 

(29) ho<C(u) Z 5 r (r -u)< H -^- a - l dudr <C{u) Z 5 r dr. 

Jo Jo Jo 

It follows from Theorem 2.2 that on Sl e ,So,p the same estimate holds for In. 
Now estimate I5. 
rt 



h<C(u) / \a(u,Y£)-a(t u ,Y°)\u- a du 
Jo 

+ C(u) f [ lair^-aiU^-v^Y^ + ait^Y^lir-uy^dudr 
Jo Jo 

— ■ Il3 + Il4- 

Obviously, 
(30) 

/13 < C(u)8P, 

ft rU ft rr 

Im< 



C(w)(y J +j J ) \...\{r -u)- a - x dudr 
<C{u) f l T 5 (r-u)- a - 1 dudr+ f f ((r - u) p + (r - u) H - p )dudr 



(31) <C(u)(5 p - a +6 H - p ~ a ). 
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Similarly, 



(32) 



h < C(u) / \a{t u ,Y°) - a(t u ,Yl)\ u~ a du 
Jo 

+ C(u) f f \a(t ri Y r s )-a(t r ,Yl)-a(t u ,Y*) + a(t u ,Yl)\ 
Jo Jo 

x (r — u)~ a ~ 1 dudr =: hs + he. 



Here 



(33) 7ib < CM / S H - p u- a du < C(uj)5 H - p ; 



(34) I 16 < C(uj) / / 6 H -»(r - u)- a ^dudr < C{uj)S H - p - a . 

Jo Jo 



Substituting (|2"T ]) -([3~4" )) into 1[2"D|). we obtain that on f2 



s,So,p 



(35) 



C(u) ( j Z 5 r r-°dr + 8 H - p - a + 5 H ~ P + j r dr) , 



where r = A r , u (X, Y s ){r - u)-"" 1 ^. Recall that H - p > 2H - 1, therefore 



CH( / (Z r 5 r-"+0. r )dr + 5 2ff - 1 -"). 



We now estimate 6* r . Evidently, for t > u 



A t , u (X,Y s )< [ \b(s,X s )-b(t s ,Y t S 3 )\ds 

J U 

m «t 

+ £ / (<T ij (^X s )-a^(t s ,Y t s e ))dBi 

i,i=i Ju 
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Therefore, using inequality ((2|), we obtain that 9t < Y2k=i ^ fc > where 
Ji= f I \b{s,X s )~b{s,Y!!)\ds{t-u)- a - 1 du, 

JO Ju 

J 2 = f f \b{8 ) Yf)-b{t t ,Y^)\da(t-u)- a - 1 du, 

JO Ju 

h= ! [ \b(t a ,Yf)-b(t ai Y t * t )\d8(t-u)- a - l du, 

Jo Ju 

J 4 = C(uj) [ [ \a(s 1 X s )-a(s,Y s s )\(s-u)- a ds(t- U y a - 1 du, 

Jo Ju 

J 5 = C(w) f f \a(8,Yf)-a(t a ,Yf)\(s-u)- a d8(Jb-u)- a - 1 d'u, 

Jo Ju 

J 6 = C{uj) f f \a{t s ,Y*)-a{t s ,Yl)\{s-u)- a ds{t-u)- a - l du, 

Jo Ju 

J 7 = CM / / / \a(r,X r )-a(r,Y r 5 )-a(v,X v ) + a(v,Y v s )\ 

JO Ju Ju 

x (r - v)- a - x dvdr{t - uj-^du, 
J 8 = C(uj) [ f f \<T(r,Y r s )-<T(t r ,Y r s )~<T(v,Y v s ) + v(tv,Y*)\ 

Jo Ju Ju 

x (r - v)- a - x dvdr{t - u)- a ~ l du, 
J 9 = C(lu) f f f \a(t r ,Y r s )-<j(t r ,Y t l)-a(t v ,Y v s ) + a(t v ,Y t s j\ 

Jo Ju Ju 

x (r - v)- a - x dvdr{t - uy^du. 

It is clear that 

Ji<C [ Z s s f (t-uy^duds, J 2 <C<5 7 , J 3 <C{u)S H - p . 
Jo Jo 

Further, 

7& 



J 4 <C Z° {s-u)- a (t-u)- a - L duds. 
Jo Jo 

As we noted before, the inner integral f*(s — u)~ a (t — u)~ a ~ 1 du < Co(t — s)~ 2a , 
= $^( l + y)~ a ~ 1 y~ a dy- Therefore 



J A <C ( (t- s)- 2a Z s s ds. 
Jo 



Similarly to J 2 , J 5 < C(uj)S~>, and similarly to J 3 , J 6 < C{uj) < C{u)5 H ~ p . Esti- 
mating J7, Jg and Jg is, of course, a bit more complicated, but not dramatically. 
Obviously, 

J%<C{u))5 p [ [ [ (r -vy^dvdrit-uy^du 



J u J u 

t 



C(w)6 p / (t - u)- 2a du < C{u>)8^; 
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similarly Jg < C(w)6 H ~ p . Now we apply to J7 the inequality (|26|) and obtain the 
following estimate of the integrand: 



(36) 



\a(r,X,,)-a(r,Y r s )-a(v.,X v ) + a(v.Y v s )\ < M |A r>l) (X, Y 5 ) 

+ \X r - Y r s \ (r - vf + \X r - Y*\ \X r - X V \ K + \X r - r/l \Y r 5 - Y v 5 \ 



According to this, we write Jj < 5~^._ 10 Jk, where, in turn, 

J 10 =C(lu) f I f A riV (X,Y s )(r - v)- a ^dv dr{t - u^^du 

J J u J u 

mv 
(t - uY a - x dul\ r ^{X, Y 5 ){r ~ v)- a ~ l dr dv 

<C(lu) f (t-r)- a 9 r dr; 
Jo 

Jh=C{lo) ( f ( \X r -Y^\{r-vf~ a ~ 1 dvdr{t-u)- a ~ 1 du 

J J u J u 

<C(lu)J^ Z s r (t-u)-"- 1 ^ (r-v) f> -"- 1 dv)dudr 

<C(lu) f (t-r)- a Z r dr, 
Jo 

J 12 =C(lu) f f f \X r -Y r s \\X r - X V \ K (r - v^^dv dr(t - u)- a - x du 

Jo J u J u 

mr 
Z s r {r - vY {H - p) - a - l dv{t - uy^dudr 
_ 

<C(lu) f Z s r (t-r)- a dr, 
Jo 

and J13 < CIuj) /„ Z*(t - r)- a dr is obtained the same way. Summing up these 
estimates, we obtain that 



J 7 <C{u) / (t-ry a (Z* + 6 r )dr, 
Jo 

whence 

(37) 9t<C(w)([ (t-r)- 2a (Z 5 r +9 r )dr + 6 H -P + <r) . 



Coupling together (|35|) and (|37|) , and taking into account that H — p > 2H — 1, 
7 > 2H — 1, we obtain 

Zf + 9 t <C(fi})($ 2H - 1 + [ ((t-r)- 2a +r- a )(Z s r +9 r )dr) 
(38) t 

< C{u)(6 2H - l +t 2a J (t- r)- 2a r~ 2a (Z s r + 9 r )dr) 

The proof now follows immediately from ([55)) and Proposition 2.1. 

The statement (ii) is obvious. □ 
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Remark 1 . In [9] it is proved that for an equation with homogeneous regular coeffi- 
cients the error \Xt — Y t s \ 5 X ~ 2H almost surely converges to some stochastic process 
£t, which means that the estimate of the rate of convergence in Theorem [5] is sharp. 

4. Approximation of quasilinear Skorohod-type equations 

Now we proceed to the problem of numerical solution of Skorohod-type equation 
driven by fractional white noise. From now on, we assume that our probability 
space is the white noise space (O,^ 7 , P) = (S'(R),B(S'(R)), fi), o is the Wick 
product, = (cj,l[ f ]) is Brownian motion, W° — B° is the white noise (see [I] 
for definitions). Next, in order to introduce an fBm with Hurst parameter H > 1/2 
on this space, we define for / : [0, T] — > R the fractional integral operator 

Mf{x)=K / {s-x) H -^ 2 f{s)ds, 

J X 

where 

K = (sin(vriJ)r(2iJ + 1)) 1/2 (Xf + K\) ~ 1/2 , 
K x = 71Y2 cos((3/4 - F/2)7r)r(3/2 - H)) , 

K 2 = 7r(2sin((3/4- H/2)tt)T(3/2 - H)) , 

and set M t (x) — Ml[ .t](x). We also define for f,g : [0,T] — > R the scalar product 
and the norm 

(f,g) H = H(2H-l) f f f(t)g(s)\t-s\ 2H - 2 dtd S , \\f\\ 2 H = (fJ)n. 
Jo Jo 

The process 

B t = (M t ,u}, te[0,T] 

is the fBm with Hurst parameter H. Let also W = B be the fractional white noise. 
Detailed description of the white noise theory can be found in [2] , [5] . 

Consider quasilinear Skorohod-type equation driven by fractional white noise 

(39) X(t)=X + [ b(s,X(s),uj)ds+ [ a(s)X(s) oW(s) ds 

Jo Jo 

with non-random initial condition Xq. Suppose that coefficients b and a satisfy the 
following conditions: 

(E) 1) The linear growth condition and Lipschitz condition on b: 

\b(t,x,w)\ < C(l + |z|), 
\b(t,x,(j) - b(t,y,u)\ <C\x-y\; 

2) "Smoothness" of b w.r.t. uj: for any t £ [0, T] and for h £ L X (R) 

\b{t,x,uj + h) -b{t,x.u)\ <C{l + \x\) I \h(s)\ds. 

Jr. 

3) Holder continuity of b w.r.t. t or order H with constant that grows 
linearly in x: 

\b(t, x, ui) — b(s, x,u))\ < C (1 + \x\) \t — s\ H ; 
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4) Holder continuity of a w.r.t. t or order H: 

W(t)-a(s)\ <C\t-s\ H . 

Remark 2. As in previous sections, we denote by C any constant which may depend 
on coefficients of the equation, on initial condition Xq and on the time horizon T, 
but is independent of anything else (and we write C(v) to emphasize the dependence 
on v). 

Remark 3. The condition (E) 2) is true if, for example, the coefficient b has stochas- 
tic derivative growing at most linearly in a;. It is obviously true if b is non-random. 

Define for t <E [0,T] <7t(s) = a(s)l[ ^(s) and denote 

J„(t) = exp° (- J a(s)dB^ = expj- J^Ma t (s)dB°(s) - i \\a t f H ds 

the fractional Wick exponent. It follows from [8, Theorem 2] that under assump- 
tions (E) equation (|39|) has the unique solution that belongs to all LP and can be 
represented in the form 

X(t) = J. lT (t)oZ(t), 
where the process Z(t) solves (ordinary) differential equation 

(40) Z(t)=X + [ J a {s)b{s,J- 1 {s)Z{s),uj + Ma s )ds. 



This gives the following idea of constructing time-discrete approximations of the 
solution of (|39|) . Take the uniform partitioning {t„ = n5, n = I, . . . , N} of the 
segment [0, T] and define first the approximations of Z in a recursive way: 



(41) 
where 



Z(0) = x , 

Z{r n+ i) = Z(r n ) + J{T n )b{r n , J~ 1 (T n )Z(T n ),uj + Ma„)S, 



J(i)=expj-^ a(s)dB s - ^ ||ffl[o,t]| 



<t(s) = a(t s ), a n = 51[o, T „]- 
Note that both ||cr ra || w and Ma n are easily computable as finite sums of elementary 
integrals. Further, we interpolate continuously by 

(42) Z{t)=X + f J{t s )b(t s J- 1 {t s )Z(t s ),LO + Ma na )ds, 



where n s — max{?i : r„ < s}, and set 

(43) X(t)=T_ Am[ot] J~ 1 (t)Z(t), 

where for h e S'(R) T h is the shift operator, T h F(u) = F(ui + h). 

Lemma 2. Under the assumption (E) 1) the following estimate is true 

\e ai b(t, e~ ai x, u) - e a2 b(t, e~ a2 x, u) | < C(l + e Ql + e° 2 + \x\) \atx - a 2 \ . 
Proof. Write 
\e ai b(t,e- ai x,u) - e a2 b(t,e- a2 x,u)\ 

< \e ai b{t,e~ ai x,uj) - e ai b{t,e" a2 x,uj)\ + \e ai b{t, e~ a2 x, uj) - e a2 b{t,e- a2 x,uj)\ 
and apply (E) 1). □ 
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Lemma 3. Let £1 and £2 be jointly Gaussian variables. Then for q > 1 



|e €l - 



2g 



<C(L, g )(E[(a-6) 2 ]) 9 , 
w/iere L = max {E [£?] , E [^1] }■ 

Proof. By Lagrange theorem, Cauchy-Schwartz inequality and Gaussian property, 
E [| e «i- e ^| 29 ] < (E[ e 4 ^ +e 4 ^] E [| Cl _ 6 |4^ 1/2 
<C(L)C(g)(E[(a-6) 2 ])\ 

as required. 

Our first result is about convergence of Z to Z. 
Theorem 3. Under conditions (E) for any p > 1 i/ie following estimate holds: 

2p 



□ 



(44) 



Z(t) - Z(t) 



Proof. Firstly, we remind that Z(t) belongs to all L q and E [ \Z(t)\ q ] < C{q). There- 
fore equation (|40| together with the condition (E) 2) gives E[|Z(t) - Z(s)\ q ] < 
C(q) \t — s\ q . Equation (|4"Tj) and the condition (E) 1) allow to write 

Z( Tn+1 )\ < (1 + C5) |z(r„)| + CSJ(r n ) < e cs \z{r n )\ + CSJ(r n ). 

This gives an estimate 

N-l 



fc=0 



Then for any q > 1 by the Jensen inequality, 

N-l 

Z(r n ) q < C(q) ]T J q (Tk)5, 



Taking expectations, we get 



fe=0 



N-l 



Z(T n )\ 9 ] < tf(«) £ E |>(t*) 



k=0 



s. 



Using that each J is exponent of Gaussian variable and a is bounded on [0, T], we 
obtain 

N-l 



Z(r n ) <C(q)J2S = C(q). 



k=0 



This through 02} and (E) 1) implies E [\Z(t)\ J < C(q). 
Now write 



Z(t) - Z(t) 



<h+h + h + h + h, 
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where 



h 



J{t s )(b(t s , J- 1 (t s )Z(t s ),u + Mo ns ) 

-b(t s , J- 1 (t s )Z{t s ),Lu + Ma ns )) ds 
(J{t s )b(t s , J- 1 (t s )Z(t s ),u J + Ma na ) 

- J a (s)b(t s , J- 1 (s)Z(t s ),co + MB n ,)) ds 



h = 



I a = 
h = 



J a {s)(b{s, J- 1 (s)Z(t s ), uj + Mo ns ) 

- b(t s , J- 1 (s)Z(t s ),u + Ma na )) ds 

f J a (s){b(s, J-\s)Z(t s ), uj + Ma n J - b(s, J-\s)Z(t s ), u + A/a,)) ds 
Jo 

f J a (s)(b(s, J- 1 {s)Z{s),^ + Ma s ) - b(s, J- 1 (s)Z(t s ), u + Ma s )) ds 
Jo 

We first estimate using Lemma [2] 

I ^ c J o {l + Ms) + J(t s ) + \Z(t s )\) (J 
+ \a(t s )(B s -B(t s )) 
<C [ (1 + J CT (s) + J(t s ) + \Z(t s )\) 



(er(lt) — cr(w)) dB v 
1 



\n _ \\ Un s \\h 



ds 



(a(u) - a(u)) dB u 



\B S - B ts \+S H )ds, 



where the inequality 



Ws\\h ~ \\VnJ\-H 



< CS H is due to E 4) and boundedness of 



a on [0,T]. Applying Cauchy-Schwartz inequality, we arrive to 



h <C 



1/2 



(1 + Jl{s) + J 2 (t s ) + Z 2 (t s ))ds 
Further, from (E) 3) 

h<C [ {J«{s) + \Z{s)\)ds5 H 
Jo 

from (E) 2) 



(B t ~B ts ) 2 +S 2H )ds 



1/2 



h <c 



f (J a (s) + \Z(s)\)ds5 H . 
Jo 



18 



YU. MISHURA AND G. SHEVCHENKO 



Condition (E) 1) allows to estimate 



h<C 



f 

Jo 



Z(t s ) - Z(t s ) 



ds, 



h<C [ \Z{s)-Z{t s )\ ds. 
Jo 

Summing up these estimates yields 



Z(t) - Z(t) 



< c 



1/2 



1/2 



+ c 



1/2 



(1 + J 2 (s) + J 2 (t s ) + Z 2 (t s ))dsj 

' 5 2H + f Q {[f Q (*(«) - *(«)) dB u) 2 + (Bt - B ts ) 2 y.s 

Z(t s ) - Z(t s ) ds + C f \Z{s) - Z(t s )\ ds. 
Jo 

Then, using (discrete) Gronwall inequality, we get 

Z(t) - Z(t)\ <c(^J o (1 + 4(a) + J 2 (t s ) + Z 2 {t s ))ds^ 

■ [s 2H + H") - *(«)) dB u) 2 + [Bt - B ts ) 

+ C [ \Z{s)-Z(t s )\ ds. 
Jo 

Then we raise this to the 2pth power and use Jensen's inequality. The last term 
will be bounded by C(p)S 2p 1 in the first one we apply Cauchy-Schwartz inequality 
for expectations, Jensen's inequality and use uniform boundedness of moments for 
Z, J fj and J (for J a and J it follows from the fact that the both are exponents of 
some Gaussian variables with bonded variance) to get 

i-T 



1/2 



ds 



Z(t) Z(t) 



Using again that E 



<C( P )(s 2 p h + (e[|jT (a(u)-a(u))dB u 



\B t - B f 



1/2 



Z(t) - Z(t) 



C(p)(E [(•) 2 ]) 2p for Gaussian variables, we get 



C(p) [s 2pH + (E [ I £ (a(u) - a(u)) dB u 

y 



2 i 



Bt~B ts \ 



< C(p)(S 2 p h + \\(t-2\\%) < C( P )5 2pH 
the last is due to (E) 4). This is the desired result. 

Now we are ready to state the main result of this section. 



□ 
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Theorem 4. Under conditions (E) approximations X defined by (|43j) converge to 
the solution X of (|39|) in the mean-square sense, and moreover 



(X(t)-X(t)f 



< cs 



2H 



Proof. Estimate first for h £ L X {R) 
T h Z{t) - Z(t) <A 1 +A 2 +A 3 



Ax 



Ao 



A, = 



T h J a {s) b(s, (T h J- 1 )T h Z(s),u + h + Ma s ) 

- b(s, {T h J- l )Z{s),u + h + Ma s ) 
ThMs) b(s, {T h J- 1 )Z{s), uj + h + Mos) 



ds, 



b(t,(T h J- 1 (s))Z(s),w+Ma s ) 



T h J a (s)b{t, (T h J~ 1 (s))Z(s), co + Ma s ) 



J a (s)b(t,J- 1 (s)Z(s),LU + Ma s 



ds, 



ds. 



The condition (E) 1) gives A x < C f* \T h Z(s) - Z{s)\ ds, the condition (E) 2) gives 
A 2 <C [ (l + \Z(s)\)ds [ \h(s)\ds 



and Lemma [2] with boundedness of a yields 

A 3 <C { (l + J a (s)+T h J(a) + \Z(s)\)ds 



Ma(s) h{s) ds 



<C (l + J a (s)+T h J(a) + \Z(s)\)ds \h(s)\ds 



R 

Applying Gronwall lemma, we get 

\T h Z(t) - Z(t)\ <C [ (1 + J a (s) + T h J(a) + \Z(s)\) ds [ \h(s)\ ds. 

Jo Jr 

Raising this inequality to the 2pth power, taking expectations and using Jensen 

inequality and boundedness of moments of Z , ,J a and Th J a (the last follows from 

the Girsanov theorem, Cauchy-Schwartz inequality and assumptions on h), we get 

\ 2p 



(T h Z(t)~Z(t)) 2p <C(p) 



\h(s)\ds 



Further, 



(X(t)-X(t)) 2 <3(A 1+ A 2 + A 3 ), 



Ax = E 
A 2 = E 
A 3 = E 



(J(t)T_ M z 1[ot] (Z(t)-Z(t))) 
((j- a (t) -J(t))T_ MS1[(Kt] Z(t)) 
(j^(t)(T_ Ma (l - T_ M ^ 1[o t] _ at) )Z(t)Y 
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where 



J{t) = cxp 



M(al m )(s)dB s ~^\\al m \\ 2 H 

Now estimate using Cauchy-Schwartz inequality, Girsanov theorem (which can be 
applied as a and a are bounded on [0, T]) and Theorem [3] 

"7*(t)]E[r_ J1 fi ?1[0it] (z(i)-z(t)) 4Vxl " 

c(E[j(t){Z(t)-Z(t)) 41 x ' 



Ax < E 



< 



< C E 



J 2 (t) E 



1/4 



< 



2ff 



Similar reasoning and Lemma [3] imply 



A 2 < CE 



J R M(al m - a t )( S )dB° s + ~( ||ot||^ - ||51 [M | 



Using condition (E) 4), we obtain A 2 < CS 2H . And for A3, using the above 
estimate, we get 



A 3 < 



This concludes the proof. 



\M(al[ ,t] ~ ct*)( s )| ds < cs 



2H 



□ 



Remark 4. It is natural to assume that the coefficient b is expressed in the terms of 
fBm B rather then in the terms of underlying Brownian motion B° (or underlying 
"Brownian" white noise u>.) This justifies the fact that it is a not Ma what is 
discretized in (|4T1). 



Remark 5. Similarly to the proof of Theorem [4] one can prove that for any s > 1 



X(t)-X(t) 



< s 



The case s = 2 is considered in the paper to keep classical "scent" of results. 

Remark 6. Results of this section can be generalized for random initial condition 
Xq in the following form: under conditions (E) and L p -integrability of the initial 
condition one has convergence in any L s for s < p with 



X(t)-X(t) 



< s 



<H 



Proofs need some simple changes: Holder inequality for appropriate powers instead 
of Cauchy-Schwartz one. 
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